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2 ' 1 Introduction 

o 
">: 

The study of nonlinear evolution equations with nondispersive travelling wave 
solutions or solitons has been at the forefront for many years now. In com- 
paratively more recent times, the field of integrable models has considerably 
S^ widened by the construction of integrable systems which are supersymmetric 

thus bringing into the supersymmetric regime many of the techniques devel- 
oped in the study of bosonic or ordinary integrable systems. Among integrable 
systems, the KdV equation [1] , undoubtedly, has been studied most intensively. 
Its supersymmetric extension in N = 1 superspace was formulated by Manin, 
Radul [2] and Mathieu [3,4] in their seminal papers. Later on, Mathieu and 
co-workers [5,6] have found integrable extensions of the KdV equation with 
N = 2 supersymmetry and more recently other supersymmetric extensions of 
the KdV have also been reported [7,8]. The integrability of these supersym- 
metric equations have been proved by showing that these equations possess 
one or more of the various criteria associated with integrability such as an infi- 
nite number of conservation laws, a bihamiltonian structure, the Lax operator 
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etc. Another aspect of integrability of a nonlinear evolution is the existence of 
soliton solutions. It has been conjectured that the existence of three soliton 
solutions implies complete integrability as it appears that the three soliton 
condition is intimately connected with the Painleve test [9]. 

A number of methods have been developed to construct soliton solutions of 
nonlinear evolution equations - the most well-known being the Inverse Scatter- 
ing Transform which was first applied by Gardner, Greene, Kruskal and Miura 
[10] to the paradigmatic KdV equation. Another technique, the Hirota bilin- 
ear method [11], is a powerful tool for computing multi-soliton solutions and 
establishing the integrability of nonlinear evolution equations. This method 
has been used by a number of workers to extract soliton solutions of many 
nonlinear dynamical systems (see for example [12,13,14,15,16,17,18,19]). The 
method involves first, in casting the nonlinear partial differential equation in 
the bilinear form by a dependent variable transformation. This procedure is 
intuitive rather than algorithmic as there is no rule for finding the correct 
transformation. As a second step, which is comparatively more straightfor- 
ward, the soliton solutions are found by substituting trial solutions for the 
new dependent variable(s) - which is (are) usually referred to as the r func- 
tion^) - in the bilinear equations and checking for consistency. Here we note 
that bilinearization of a nonlinear evolution implies the existence of one and 
two soliton solutions but three and higher solitons may not exist. It is the 
three soliton solution whose existence is crucial for integrability. 

The Hirota formalism has also been adapted to bilinearize supersymmetric sys- 
tems [20,21,22,23,24,25] and soliton solutions of a number of supersymmetric 
integrable equations have been obtained by the application of this technique. 
For instance, the N — 1 KdV of Manin-Radul-Mathieu [2,3,4], 

d t 4> + ^D 2 (4>D(f)) + D 6 <j> = (1) 



was bilinearized [20] via the dependent variable transformation 

= 2D 3 logr (2) 

In (1) and (2), D = dg + 9d x is the supersymmetric derivative, and 9, the 
Grassmann odd variable. The transformation (2) allows one to write the N — 1 
KdV equation in the bilinear form [20] 

(SD t + SD3)(r.r) = (3) 

where S is the supersymmetric generalization of the Hirota operator in N = 1 



superspace and is denned by [20] 



SB n x (f.g) = (D 01 - D e2 )(d xi - d X2 ) n f(x 1 ,6 1 )g(x 2 ,6 2 ) 



(4) 



X\ — x 2 — x. 
0i = 2 = 



and D is the ordinary or bosonic Hirota operator given by 
D "(/-^) = (d Xl -d X2 ) n f{x 1 )g{x 2 )\ Xl=X2=x . 



(5) 



The soliton solutions for the N — 1 KdV were found by the Hirota method 
[21,22] through the standard procedure of substituting trial functions for the 
r-functions and then checking for consistency. It was observed that the two and 
higher order soliton solutions could be found only if constraints are imposed on 
the fermionic parameters. Subsequently, Carstea, Ramani and Grammaticos 
[24] have shown that two and three soliton solutions for the N — 1 KdV can 
be obtained without imposing additional fermionic constraints by dressing the 
fermionic part of the soliton solutions through the two soliton interactions and 
these solutions, therefore, constitute a more general class of soliton solutions. 

In this paper, we have attempted to find if there are other supersymmetric 
nonlinear evolution equations which have general solutions of the type con- 
structed in [24]. Earlier work on the bilinearization [23] of the N — 1 mKdV 
[3,4] provides the clue to finding one such system. The N = 1 mKdV 



dtip + D & 4) - 3ijD 3 i;Di) - 3(Dif)) 2 D 2 ij = 



(6) 



can be converted into the following set of four bilinear equations [23] 
(SD t + SD|)(n.r 1 ) = 
(SD t + SD3)(r 2 .r 2 ) = 
SDJn.ra) = 



(7) 



(9) 



and 



D'(r!.r 2 ) 



(10) 



after transforming the superfield ip to 



i) = D lo§ 



:n) 



r 2 



It was noted in [23] that for the N — 1 mKdV, in one of the bilinear forms, 
namely (10), the supersymmetric Hirota operator S does not appear. As a 
result, two and three soliton solutions without constraints on the fermionic 
parameters could not be obtained for it. However, if a supersymmetric system 
which had bilinear forms (7), (8) and (9) only could be constructed, one would 
have constraint-free soliton solutions of the type obtained in [24] . 

In the next section, we construct a supersymmetric system in N = 1 su- 
perspace that can be transformed into the bilinear forms (7), (8) and (9). 
Moreover, the supersymmetric system has the KdV as its bosonic limit and is 
therefore a supersymmetric extension of it. In Section 3, one, two and three 
soliton solutions for it are written down. These soliton solutions are of the type 
found in [24] i.e. free from constraints on the fermionic parameters. Section 4 
is the concluding one. 



2 An N = 1 Extension of the KdV 

Consider the following pair of supersymmetric equations in N = 1 superspace: 
d t cj> + dl<P - 6d x [(D(f>)<f> £>"V] - 2d x cf> 3 + 3^(^00) = (12) 

dtifj + d\i> - Qd x [(Di/j)<p D"V] - Zd^D^) + Zd x [{dx<l>)D<i>) = (13) 

where D~ l = 9 + dgd~ l is the formal inverse of the operator D i.e. DD^ 1 = 
D~ l D = 1. Clearly, (12,13) are nonlocal, coupled, partial differential equations 
in N = 1 superspace for the superfields <p and ip which have conformal weights 
1 and | respectively. If these superfields are written in terms of their bosonic 
and fermionic components 

4> = (j) b + % (14) 

il> = ij f + e^ b (is) 



then (12) and (13) can be split up into four equations in terms of the compo- 
nent fields. For (12), they are 



dtfo + dl<p b - 64(0 /( /> 6 cTV/) - 2d x $ + 33,(^/0/) = (16) 

dt<f> f + d 3 x <p f - 69 x [(9 x ^)0 6 9jV/] + ZdxiMf) - Zd^fdM = (17) 

and the component equations for (13) are 



dtipf + &j>s - 6^(^0 6 a;V/) - zd x ($ffo) + w x {(t>sd x <j> h ) = o (is) 

d t ip b + <%ij) b - Qd x (ip b 4> 2 b ) - 6ip b d x i/j b + 3d x (i/j f d x ip f ) + 3d x (d x (p b ) 2 
-Gd^^fd- 1 ^) - 6^[(^/)<^-V/] " 34(0/^0/) = (19) 

Equations (16), (17), (18) and (19) are invariant under the N — 1 supersym- 
metric transformations 



5(fr b = rjcpf 5<p f = rjd x (j) b 
5iP b = r]d x i) f 8ip f = rii> b 



(20) 
(21) 



and therefore constitute an N = 1 supersymmetric system. Moreover, on 
setting the fields <p b , 0/ and ipf to zero we obtain the KdV equation 



d t tp b + 8%i/) b - 6ip b d x ip b = 



(22) 



This clearly indicates that equations (12,13) may be considered a supersym- 
metric extension of the KdV equation. Note that, as in case of the N = 2 KdV 
[5], the bosonic core of the equations (12) and (13) also contain the mKdV 
equation d t (f) b + <9 3 b — 2d x §\ = 0- 

In order to write the supersymmetric equations (12) and (13) in bilinear form, 
consider the following transformations for the superfields and i[): 



D 2 log 



n 

r 2 



(23) 



and 



^ = J D 3 log(r 1 r 2 ) 



(24) 



where T\ and r 2 are independent, bosonic superfields. Substituting for the fields 
and ip in the equations (12) and (13), we find that they can be cast into 
a form involving only bilinears. For instance, equation (12) after substitution 
acquires the form 



d t D 2 log — + D 8 log — - 6D 2 

r 2 r 2 

3 



D 3 log— £> 2 log— Dlog — 



72 









-2D 2 



D z \og 



T 2 



+ 3D' 



L> d log(nr 2 )Dnog 



T 2 



(25) 



Adding and subtracting the term 6D 4 log— D 3 log(r 1 r 2 ) to the left hand side 
of (25) after operating on it by D^ 1 we get, 



2d t D log n + 2D 7 log 7i + 12L> 4 log n D 3 log n 
2a t D log r 2 + 2D 7 log r 2 + 12D 4 log r 2 D 3 log r 2 



-12(£> 4 logri-D 4 logr 2 



n 



T\ 



£> J log(riT 2 ) + D 2 log-Dlog 

T2 r 2 



(26) 



Writing the above equation using the Hirota operators, we find that equation 
(12) is finally transformed into 



(SD, + SD3)(n.ri) (SD t + SDf)(r 2 .r 2 ) 



rf 



D^(ti.Ti) D2(r 2 .r 2 ) 



T 2 
r 2 



?2 

SD a (ri.r 2 ) 



(27) 



An almost identical calculation transforms equation (13) into 



(SD t + SD3)(n.ri) , (SD t + SD3)(r 2 .r 2 ) 



- 6 



SD x (r 1 .r 2 ) 

T\T 2 






D^(r 2 .r 2 ) 







(2S 



From (27) and (28), it follows that the supersymmetric system (12,13) has the 
bilinear forms given by (7), (8) and (9), that is only those equations which 
involve the super Hirota operator S. Clearly, unlike the N = 1 mKdV equation, 
this supersymmetric system does not have any bilinear equation involving only 
bosonic Hirota operators. 



3 Soliton Solutions 



For one soliton solutions, consider the following forms for the r functions: 



T\ 



D v+ec 



(29) 



and 



r-2 



?V+8( 



(30) 



where rj = kx + out. Here k and ui are the bosonic parameters while £ is 
the Grassmann odd parameter. These are solutions of (7) and (8) provided 
u = — k 3 which is the KdV dispersion relation and in addition, the r functions 
(29,30) satisfy (9) and are therefore appropriate r functions for the set bilinear 



equations (7,8,9) representing (12,13). Throughout this section, we employ the 
notation similar to that used in [24] for the r functions of the N — 1 KdV. 

In obtaining the two soliton solution for the N = 1 mKdV equation [23], it 
was found that the solutions exist provided the fermionic parameters £i and 
(2 are related via the constraint relation k\C, 2 — &2C1 an d identical relations 
characterize the three soliton solutions also. For the supersymmetric equations 
(12,13) however, by an appropriate choice of r functions, general solutions can 
be obtained. Explicitly, the r functions are 



T\ 



P m+0Ci 



1 + e" 11 ^ 1 + e rm ' 2 + A 12 



,m+0(2 



1 



C1C2 



k 2 - k\ 



D »?l+»72+0(ai2<l+a2lC2) 



(31) 



and 



r 2 



3 ^i+eci 



D m+sc,2 



A 



12 



C1C2 



k 2 - h 



ril+ r l2+8(ai2(i+a2i(2) 



(32) 



where 77* = kiX + corf, (i — 1, 2) and 

2 



A 



12 



h - k 2 
h + k 2 



(33) 



is the bosonic interaction term. The a^ 's which are responsible for the dress- 
ing of the fermionic part through the two soliton interactions [24] are given 
by 



a.-. 



ki T" Kj 



(«,J' = 1,2). 



(34) 



Substituting (31) in (7) and (32) in (8) yields the dispersion relations Ui = 
—kf with % — 1,2. These r functions also satisfy (9) and are therefore the 
appropriate r functions for the system considered here. 

The r functions for the three soliton solutions of the system (12,13) have the 
form 



Tl = 1 + e Vl+eCl + e m+eC2 + e r?3+9C3 



^l+»?2+6(ai2Cl+ a 2lC2) 



+A 12 


[l + 2, ClC2 , 1 

«2 - «1_ 


+A n 


[l + 2/ ClC3 , 1 
«3 - «4_ 



a «7lH-»/3+0(ai3Ci+Q:3lC3) 



and 



+A 



23 



1 + 2 



C2C3 



fca - fe 



3 >72 +»?3 +^(«23 C2 +«32 C3 ) 



+A 12 A 13 A 



12^ 1 13^ 1 23 



1 + 2 



(ttl3Cl)("23C2) 



fc 2 - fci 



1 + 2. 



^i 2 Cl)(«32C3 
fc 3 - fei 



1 + 2 



^2lC2)(Q3lC3 



»)l+»?2+r;3+^(oi2Ol3Cl+ a 21O23C2+a3ia32C3) 



ri = 1 — e 



^i+^Ci _ ^2+0<2 _ „r? 3 +6»C3 



— e 



+^12 
+A 13 

+A 23 



1 + 2 
1 + 2 
1 + 2 



C1C2 



-AwAkA 



12^113^-23 



k 2 - h 

C1C3 
h ~ h 

C2C3 
k 3 - h 

1 + 2 



o m+V2+9(otl2(l+a-2l(2) 



rii+ri 3 +6(a 13 (i+a 31 ( 3 ) 



3V2+m+8( a 23(2+OC32C3) 



(ai 3 Cl)(«23C2 



&2 - h 



1 + 2 



>l 2 Cl)(«32C3) 

fc 3 - k t 



1 + 2 



(ogiCgKogiCg) 

^3 - k 2 



o rjl+V2+m+0( a 12Otl3Cl+a21Ol23C2+a3ia32C3) 



Here the bosonic interaction terms A- are given by 



(35) 



(36) 



A 



1,1 



, ^i + fy , 



(i,j = 1,2,3) 



(37) 



and the a^ 's are obtained from (34) with i,j = 1, 2, 3. We have verified that 
these r functions satisfy the bilinear equations (7), (8) and (9) and hence 
three soliton solutions exist for the supersymmetric system given by (12,13). 
Also note that function T\ for the one, two and three soliton discussed above 
are exactly those found by Carstea-Ramani-Grammaticos [24] while the cor- 
responding t 2 's differ by the sign of some of the terms. The relation between 
structures of the r functions for the N — 1 KdV [24] and those written down 
here parallels the relations between the forms of the r functions for the bosonic 
KdV [11] andthemKdV [12]. 

Thus the supersymmetric evolution equations (12,13) have at least three soli- 
ton solutions which are free from fermionic constraints. Since the existence 
of three soliton solution implies integrability, (12,13) constitute an N = 1 
supersymmetric integrable extension of the KdV equation. 



It is also relevant to mention that for extensions of the KdV with a single 
additional fermionic field, there are only two which are know to be integrable - 
the super KdV of M anin-Radul-Mathieu [2,3,4] which is invariant under N — 1 
super symmetric transformations and the one proposed by Kupershmidt [26] 
which is not. Among N = 2 extensions, i.e. those with two fermionic and an 
additional bosonic field and also invariant under N = 2 supersymmetry, five 
are known to be integrable. They are the three N = 2 super KdV equations 
which are integrable for a free parameter a having values 2, —1, 4 [5,6], the odd 
N = 2 super KdV [7] and the super KdV-B [8]. Although system consisting 
of equations (12,13) discussed in this paper is a superextension of the KdV 
with two fermions and an extra bosonic field, it is however not invariant under 
N = 2 supersymmetry but is an N = 1 supersymmetric, integrable extension 
of the KdV. A distinctive feature of this system is its nonlocality which is not 
observed in the integrable extensions of the KdV that are already classified. 



4 Conclusion 



In this paper, we have constructed a pair of coupled N = 1 supersymmet- 
ric evolution equations whose bosonic limit contains the KdV equation. This 
system was constructed in such a way that it can be transformed into a set 
of Hirota bilinear equations which possess at least three soliton solutions free 
from constraints. Since three soliton solutions exist for the system, the system 
may be considered an integrable, N — 1 supersymmetric extension of the KdV 
equation. 
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